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Abstract. Crystal plasticity is the result of the motion and interaction of dislocations. There is, however, still a major 
gap between microscopic and mesoscopic simulations and continuum crystal plasticity models. Only recently a higher 
dimensional dislocation density tensor was defined which overcomes some drawbacks of earlier dislocation density measures. 
The evolution equation for this tensor can be considered as a continuum version of dislocation dynamics. We use this 
evolution equation to develop evolution equations for the total dislocation density and an average curvature which together 
govern a faithful representation of the dislocation kinematics without having to use extra dimensions. 
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INTRODUCTION 

Dislocation based continuum modelling of plasticity started in the 1950s with the works of Nye [? ], Kroner [? ], 
Kondo [? ] and Bilby et al. [? ], who introduced largely equivalent tensorial descriptions of the dislocation state of a 
crystal. It was soon recognised that these concepts were not suited as foundations of a plasticity theory. This, together 
with their reasonably high mathematical prerequisites, condemned these approaches to a rather shadowy existence 
since ever. They were revived in the 1990s after Fleck, Ashby et al. [? ] recognised the role of the Kroner tensor in 
(size-dependent) strain gradient theories as a measure of so-called geometrically necessary dislocations [? ] (GND). 
Since then various strain gradient dependent plasticity theories have been developed, which are sometimes combined 
with evolution equations for the so-called statistically stored dislocations (e.g. [? ] [? ]) (SSD). The latter evolution 
equations are usually based on a second direction of dislocation based modelling which traces back to Kocks [? ]. 
He was the first to introduce phenomenological evolution laws for the (scalar) total dislocation density. However, 
the available evolution equations for SSDs are usually driven by plastic strain and local in the sense that they do 
not account for dislocation fluxes. This may be justified for polycrystals where dislocations mostly accumulate strain 
driven because they are trapped at grain boundaries. On the single crystal level, however, plastic strain accumulates 
where dislocations have passed through, while the dislocation density accumulates where dislocations get trapped, i.e. 
where plastic strain is low. 

In [? ] we showed that dislocation fluxes and the evolution of plastic strain can be handled consistently by using 
a higher dimensional dislocation density measure. In the current work we show how the extra dimension may be 
eliminated from the evolution equations by integration to yield evolution equations for the total dislocation density 
and the average curvature. These evolution equations are able to handle dislocation fluxes and line length changes 
without using extra dimensions. 

EXTENDED CONTINUUM THEORY OF DISLOCATIONS 

At the heart of the extended continuum theory of dislocations developed by the authors [? ] lies the so called dislocation 
density tensor of second order (SODT) a 11 . This tensor is a natural generalisation of the classical dislocation density 



tensor to a higher dimensional configuration space. It is closely related to the phase space densities of dislocations as 
introduced by El-Azab [? ]. 

If dislocations move by glide only, the SODT is defined on the configuration space Q = M x \Jp Sp , where M denotes 
the (spatial) crystal manifold, j3 indicates the slip systems and the are unit circles of directions in the respective 
glide planes. A point in Q is considered as composed of a spatial point p and the angle <pa between a direction in the 
glide plane and the Burgers vector bp . We may consider the SODT as a sum of tensors defined for each slip system, 
i.e. a 11 = CXp- On each slip system, the SODT is defined by a density function pp (p,(pp) giving the average number 
(per unit area) of dislocations at p with line direction lp((pp), and a curvature function kp(p.<pp) characterising the 
average curvature of these dislocations. In the following we restrict ourselves to a single slip system (and skip the index 
/3) and chose a coordinate system such that the slip plane is the 1-2 plane and the Burgers vector points in 1 -direction: 
b = b x d\. The canonical line direction consequently reads l((p) = (cos <p,sin<p). The second order dislocation density 
tensor then takes the form 

a 11 = p (p,(p) (cos (p, sin (p,k(p,(p)) ®b. (1) 

We call L = (cos (p, sin (p,k(p,(p)) the generalised line direction. Because dislocations do not end inside a crystal the 
SODT needs to be solenoidal, that is 

Diva 11 = (cos<pdip(p,<p) + sin(pdip(p,<p) + d <p (pk(p,<p)))®b = 0, (2) 

where Div denotes the divergence operator on the configuration space. We note that the last equation is equivalent to 
Div(pL) = 0. Note that in the sequel we will usually drop the arguments of the involved objects to maintain readability. 

We assume that the dislocations move in the direction of the spatial velocity field v = (v',v 2 ) :— 
v (p,(p) (sin (p,~- cos<p) perpendicular to the line direction I. Furthermore we introduce the generalised velocity 
V = (v 1 ,v 2 , V_lv), where the third component Vj,v reflects the rotational velocity of a moving dislocation segment. 
The evolution equation of the SODT may be given as evolution equations for p and k (cf. [? ]) which read 

d t p = -Div(pV)+pv = -div(pv) + <9, p (pV z v) +pvfc and (3) 
d,k = -vk 2 ~V L V L v + V v k. (4) 

For later use we note that instead of the evolution equation for k we may also look at the evolution of the product pk. 
The evolution equation for pk is easily obtained from (01 and as 

d t (pk) = -Div (pkV)- Div (V £ vpL), (5) 

under consideration of the solenoidality of a 11 , i.e. of Div (pL) = 0. 

Relation to the classical dislocation density measures 

The total dislocation density p t can be determined from a 11 as 

(■In 

p t = / p((p)d<p. (6) 
Jo 

The classical disclocation density tensor a is obtained by integrating the density against the line direction, i.e. through 

a= p ((p) (cos (p, sin (p)d(p®b. (7) 
Jo 

As we consider only one Burgers vector we can decompose the dislocation density tensor into a = K®b with a vector 



k=(k 1 1 k 2 ):= p (<p) (cos<p,sin<p)c/<p. 
Jo 



(8) 



CONTINUUM THEORY OF DISLOCATIONS WITHOUT EXTRA DIMENSIONS 



Under the assumption that the (scalar) velocity v does not depend on the dislocation line direction, the evolution of the 
total dislocation density is easily obtained from Eqn. (f3]) as 

d t pt = J (-dW(pv)-d (p {pV L (v)) + pvk)dq) (9) 
= -div J p (sin(p,- cos (p)d(p^j + vj pkd(p (10) 

= -&iv{v{K 2 ,-K l ))+V J pkd(p. (11) 

The first term on the right hand side measures the net dislocation flux. That this is the right flux term to describe the 
change in total dislocation density can be seen by realising that dislocations with opposite line directions also move in 
opposite directions. Therefore only the difference between the amount of dislocations of each 'sign' leads to a change 
of total line length in a volume. The second term accounts for the change of total line length connected to the expansion 
or shrinkage of curved lines. 

As in the extended continuum theory the evolution equation for the total dislocation density needs to be accompanied 
by an evolution equation for the (average) curvature k in order to capture the line length changes correctly. We make 
the following definitions: 

pk := / pkdcp (12) 
Jo 

k := Pl^lilPMl (13 ) 
Pt S^pdcp 

With these definitions and the additional notation K 1 - — (ff 2 , — K l ) we rewrite Eqn. (fTTTt as 

d,p t = -div fvK" 1 ) + vpk = -div (vK^J + vp t k. (14) 

To derive a closed set of evolution equations we first look at the evolution equation for the averaged product pk. 
Note that as in Eqn. (O the directional part of the generalised divergence Div vanishes upon integration over <p and we 
find 

rill 

d t pk = - (div(pfcv)+div(V L v-p/)W<p (15) 
Jo 

= -div J p k(sin q>,- cos (p)d(p^J -div fj V L v ■ p (cos q>, sin (p)d(p^J . (16) 

In the evolution equation for p t the divergence term could be expressed in known quantities under the assumption that 
the velocity v does not depend on the line direction. In Eqn. (TT3T > we can close the first divergence term by assuming 
that also the curvature is the same for dislocations of each direction. This is a reasonable assumption in quasi-static 
situation where the curvature essentially balances the local shear stress. With this assumption we find for the first 
divergence term 

div J pk(sin(p, - cos (p)d(p^j — div ^vkK 1 -^ . (17) 

for the second divergence term in Eqn. ( 031 ) we need a further simplifying assumption. Here we assume that the 
dislocation density is nearly independent of the orientation which implies that the amount of GND is small compared 
to the total dislocation density, ||fc|| << p t . Note, however, that a small amount of GNDs does not imply an equal 
distribution of density in the orientation space. We think that a nearly isotropic orientation distribution is a common 
situation when dealing with mesoscopic average volumes. In this case p is largely independent of <p and we necessarily 
have p(<p) = Pt/27T. In the following we assume that the latter holds as equality. Furthermore we again assume v to be 



isotropic and find for the second divergence term 



<& v (j ^lv- p (cos (p, sin (p)d<p\ = div J (cos(pdiV + shi(pd2v)(cos<p,sm(p)d(p\ (18) 

= div(y t (d lV ,d 2 v)), (19) 



V 

where we used 

fin fin rln 



cos ((p)d(p = / sin ((p)d(p = % and / cos (<p) sin <p<i<p = 0. (20) 
Jo ' Jo 

In the case of v, p and k each being (nearly) isotropic we consequently have 

d t pk = -div (v~kK^ - ^div (p,Vv) . (21) 

As in the case of the extended continuum theory we may alternatively look at the evolution of k which we can now 
easily determine from its definition and basic differential calculus as 

d,k = d t (^] (22) 



= _vF-I(vV r x*+iv V vPtJ-^Av. (23) 

Here A denotes the two-dimensional Laplace operator. 

The evolution equation for the dislocation density tensor was derived in [? ] and translates to the evolution of K as 

d,K - (dziptvl-diiptv)). (24) 

Under the premise that the assumption of isotropic p, v and k holds, Eqns. ( fT4l . j23l , and (1241 thus define a 
kinematically closed system of equations for the evolution of an averaged dislocation system. Note that a fully closed 
system would need a relation between v, the applied stress and the current dislocation state. 



SUMMARY AND OUTLOOK 

In the current work we used a higher dimensional continuum theory of dislocations to derive coupled evolution 
equations for the total dislocation density, the classical dislocation density tensor and the average curvature of 
dislocations. The most important novelty in these equations is that they take into account the change of total dislocation 
density due to dislocation fluxes. To describe the evolution of the average curvature several simplifying assumptions 
are needed which we seek to weaken in future work. 



